Driven diffusive system: A study on large n limit 
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We study the generalized n component model of a driven diffusive system with annealed random 
drive in the large n limit. This non-equilibrium model also describes conserved order parameter 
dynamics of an equilibrium model of ferromagnets with dipolar interaction. In this limit, at zero 
temperature a saddle point approximation becomes exact. The length scale in the direction trans- 
verse to the driving field acquires an additional logarithmic correction in this limit. 
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I. INTRODUCTION 

In nature, a large class of phenomena are far-from- 
equilibrium in origin. Unlike the equilibrium systems, de- 
pending on the complexity and the relevant questions in- 
volved, there are different ways to treat non-equilibrium 
systems. In equilibrium, due to the random energy im- 
parted by the environment (thermal noise), all possible 
energy states are accessible and there is no net flow of 
probability (i.e. no net current) from one configuration 
to the other. As a result there is a straight forward way of 
obtaining the free energy and hence the thermodynamic 
properties by just summing over all the Boltzmann fac- 
tors exp[— E/KbT] for all the states, where E is the en- 
ergy of the state. In general, although one may write 
down a master equation describing the time evolution of 
the probability of a system being in a state, in the large 
time limit, the equilibrium dynamics ensures a Boltz- 
mann distribution. This behavior is preserved by the so 
called detailed balance principle. This shows up also in 
the Langevin equation which describes the time evolution 
of the degrees of freedom (say e.g. the order parameter in 
a spin system, or the height of a growing surface). In the 
equilibrium case, the Langevin equation can be obtained 
from the Hamiltonian representing all the relevant inter- 
actions of the system. However, in the non-equilibrium 
cases, such a description through a Hamiltonian is not 
possible due to one or more terms. These terms are 
purely of non-equilibrium origin and appear due to, say, a 
finite driving force in the driven dynamics of particles Jl]] 
or in random deposition of particles in the case of crys- 
tal growth ||. There are non-equilibrium systems where 
detailed balance may also hold good. In these cases the 
dynamics reduces to the equilibrium dynamics of a sys- 
tem described with appropriate Hamiltonian [^) . 

These aspects of non-equilibrium dynamics are present 
in different variants of driven diffusive systems (DDS) 
where particles with diffusive dynamics have biased hop- 
ping due to a driving force in a particular direction. In 
the case of fixed time independent driving force, a formu- 
lation of the Langevin equation through a Hamiltonian 
is not possible due to the term responsible for the drive. 



This is not necessarily true for all the driven systems. 
In particular, we shall later be concerned with a driven 
system with an annealed random drive where the prob- 
lem effectively reduces to the equilibrium dynamics of a 
specific magnetic system. In general, without the driving 
force, the equilibrium system is like a lattice gas where 
particles randomly occupy lattice sites and also due to a 
nearest neighbor attraction a particle prefers to have its 
next neighboring site occupied by another particle. In 
this sense the equilibrium static properties are the same 
as that of the Ising model with a ferromagnetic interac- 
tion. Even though, the dynamics, unlike the spin case, 
is different due to the total particle number conservation 
in the lattice gas, the static properties of the two mod- 
els are the same due to the equivalence of canonical and 
grand-canonical ensembles. With a finite driving force 
||, the dynamics of lattice gas is significantly different 
from the equilibrium order parameter conserving dynam- 
ics. Three key features, such as, the boundary condition, 
particle number conservation and the driving force to- 
gether affect the dynamics in a nontrivial way and as 
a result one observes generic power law correlation far 
above the criticality and anisotropy in the exponents of 
different correlation functions Q. 

Due to the time dependent probability distribution, 
and lack of detailed balance, many questions remain un- 
resolved for this system. Although in the steady state, 
where the distribution settles to a time independent form, 
one can write down a local detailed balance form for 
the transition rates in the Master equation, obtaining 
the physical properties from this is, in practice, difficult. 
However, even working with very small lattices, certain 
interesting features can be observed For example, 

on a small size half filled lattice, in the ground state, un- 
like the equilibrium case, one can see a strip-like occupied 
region in the presence of a large driving force. A care- 
ful analysis of the internal energy fluctuation and specific 
heat also reveals new consequences of the driving force 

In view of the complications, it is, therefore, often 
useful to adopt different approximation schemes or nu- 
merical techniques. Since in many cases the interest is 
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mostly on long wavelength and long time features, the 
continuum formulation in terms of the Langevin equation 
comes quite handy. The continuum formulation provides 
a fruitful basis for performing the mean field type calcu- 
lations or field theoretic renormalization group analysis 
|| with the latter being especially successful in explain- 
ing the new features at the critical point. In this scheme, 
the existence of a non-equilibrium stable fixed point that 
governs the new universal properties, anisotropic non- 
Gaussian scaling exponents associated with it, is clearly 
observed. The fact that far above the criticality, the two 
point or three point correlation exhibits a power law be- 
havior, can be directly shown using the above continuum 
formulation and averaging with exp[— J\ with J as the 
dynamic functional. This weight factor essentially plays 
a role that is analogous to the Boltzmann weight for the 
equilibrium system. 

Although significant progress has been made in under- 
standing the properties far above and at criticality, the 
situation below criticality is relatively unclear. Below the 
critical temperature, one would expect phase separated 
occupied and unoccupied regions provided the filling is 
not too far away from the half filling where one would 
remain out of the coexistence regime. In the presence 
of the driving force, the shape of the coexistence curve 
changes from the equilibrium situation and also the 
shape of the phase separated regions gets distorted. The 
basic questions that still remain are related to these as- 
pects. Although numerical simulations show that the 
particle occupied regions are somewhat elongated in the 
direction of the field, sufficient analytical progress is yet 
to be made. 

Along with this development, several variants of this 
simple case of driven lattice gas have appeared in order 
to understand more realistic phenomena in nature. One 
simple choice of the boundary condition is the periodic 
one in both transverse and longitudinal direction. In that 
case the system looks like a torus with the driving field, 
possibly a uniform electric field looping around it. To ob- 
tain such a situation in laboratory, one would require a 
magnetic flux increasing linearly in time. Therefore due 
to these practical purposes, a demand for more abundant 
situations with non-uniform drive or different boundary 
conditions was obvious. Thus came up a simpler but 
more realistic model of DDS with annealed random drive 
Interestingly, this model exhibits new universal 
properties that are different from the uniformly driven 
Ising model. 

In this paper, we arc interested in studying a driven 
diffusive system with annealed random drive with n com- 
ponent order parameter in the large n limit. In this limit 
the problem is exactly solvable at zero temperature. It is 
known that the large time behavior of a system, quenched 
below the critical temperature, is characterized by the 
equal time structure factor 



c(M) - (&(M)&(-M)> = L a (t)F(kL(t)), (l) 

where L(t) ~ t 1 ^, with z being the dynamic exponent, 
represents the characteristic size of the correlated regions 
growing with time fiofl . In the above definition, the sub- 
script i denotes the component of the order parameter 
field 4>. This scaling law is preserved in the case of 
quenching of a system with a non-conserved order pa- 
rameter. However, this scaling form is not maintained 
in a system with a conserved order parameter. In the 
latter case one has two length scales, obeying two differ- 
ent scalings Jll],[l2j ■ While the growth of one length scale 
with time is governed by the usual dynamic exponent, 
for the other length scale, an additional time dependence 
appears. The reason for such multi-scaling is understood 
to be the non-commutativity of the large n and large t 
limits H). 

II. MODEL 

In the lattice gas systems with a fixed total number 
of particles, there is an underlying continuity equation 
obeyed by the particle density p, d t p = —V • J, where J 
is the current. In the absence of the driving force, the 
current can be expressed as the gradient of the chemi- 
cal potential — AV/x, where A is the transport coefficient, 
p = SH/Sp, with 7i as the appropriate lattice gas Hamil- 
tonian. With an attractive interaction among the par- 
ticles, the complete description of the dynamics involves 
only one field, the particle density which plays the role 
of the order parameter. Near equilibrium, the equation 
of motion for the conserved dynamics of such attracting 
particles can be written using the continuity equation and 
the Hamiltonian, expressed in terms of the local magne- 
tization 0(x, t) — 2p(x, t) — 1, as 

H{<f>) = J d d x{\{Wct>f + T -<? + -^ 4 }- (2) 

The drive, which is actually responsible for the far-from- 
equilibrium behavior, gives rise to an additional cur- 
rent, which must vanish if there is no hole or no par- 
ticle. Thus, the simplest choice for the current may be 
J e = 4p(l — p)E = (1 — 4> 2 )E, where E is the strength 
of the driving field. The system is anisotropic due to the 
driving force in a particular direction, henceforth denoted 
by z. In the equation of motion, this z direction and the 
transverse d — 1 dimensional space are distinguished by 
different coefficients as shown below 

dttfat) = \{{t ± - V 2 )V 2 + (t,| -a\\d 2 )d 2 ^ 
-a x 9 2 V 2 + ^(V 2 ^ 3 + X d 2 4?) + Edcf, 2 } - (V • £ + dQ, (3) 

Here and in the following discussion, V in general ex- 
tends over the d — 1 dimensional space and d denotes 
the derivative with respect to z only. The last two terms 
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in the Langevin equation originate from the usual noisy 
part of the current and have short range correlations both 
in space and time with variances n± and nn. The co- 
efficients associated with the derivatives with respect to 
transverse or parallel coordinates have subscripts || and _L 
respectively. As the critical temperature is approached, 
the most realistic choice of parameters is tii > and 
t± — > [Q. In fact, it can be explicitly shown that while 
the transverse direction, and hence t± , remain unaffected 
by the driving force, rii is renormalized upwards due to 
the driving force and remains uniform for the whole sys- 
tem. 

The focus of the present paper is on an annealed ran- 
dom drive which has a correlation (E(x.,t)E(x.' ,t')) — 
<7<5(x — x.')S(t — t 1 ). A well-known starting point for the 
dynamics is to write down the dynamic functional J that 
involves the Martin-Siggia-Rose response field <f> [[l5) . Af- 
ter averaging out the drive, and ignoring all the redun- 
dant or irrelevant terms, the generating functional can 
be written as 



Jb 



dtd d xX{4>(x,t)[\- 1 d t - (r - V 2 )V 2 - <9 2 ]</>(x, t) 
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+ <t>V<f>}. 



(4) 



Since the non-Gaussian interaction term has no longitu- 
dinal operator, it does not affect rii, which, as a result, 
has been considered to be a constant (unity) in (Q) . The 
coefficient of the last term in the curly bracket has also 
been set to unity by suitable reparametrization. For con- 
venience we have denoted tj_ in this equation by r. 

It is now appropriate to point out that there exists a 
direct connection between the driven system and uniaxial 
ferromagnets with dipolar interaction [^],[l6 17 1. The gen- 
erating functional can be written in the detailed balance 
form as 



J b = / {4>d t 4>- 4>w 



where 



hm = y fe -</>(-fc)5 o - i (fc)0(fc) + 

4>(k 1 )4>(k 2 )(l)(k 3 )(j>(-k 1 -k 2 - k 3 ), 



(5) 



(0) 



hi ,/C2 )&3 



with Sq (k) = fcj^rfc 2 + k\ + fc|]. This Hamilto- 
nian represents the static system of uniaxial ferromag- 
net with dipolar interaction Jl^Jl^]. Therefore studying 
this driven system with annealed random drive also im- 
plies studying the order parameter conserved dynamics 
of uniaxial ferromagnets. 



III. RESULTS 

Our starting point for the large n analysis is essen- 
tially an n component generalization of the Hamiltonian 



Tid appearing in the detailed balance form above. To pro- 
ceed further, we start with the Langevin equation which 
represents the dynamics of the driven diffusive system 
of our interest as well as the conserved dynamics of the 
uniaxial dipolar system 



at 



- =V Q [V a |5+C«(M)], 



(7) 



where ( a (x, t) represents the short range correlated Gaus- 
sian noise. Here a (— 1, 2, ...d) stands for the dimension 
of the space and i (— 1, ....n) denotes the order param- 
eter component. For the present discussion, we shall re- 
strict ourselves to the zero temperature situation where 
we need not bother about the noise term. 

In Fourier space the, Langevin equation becomes 

^ = {-rkl -ki- k\)^ - |a(i)fti^, (8) 



where 4> 2 = Y%=i 



obta 



f) = o(t). Solving (§, we 



<f>i(k,t) = l (A: ; O)e- (T ' £ i + ^ + ' £ ii )t - (tlo/3!) ^ h(t) , (9) 
with&(t) = f*a(t')dt'. Therefore 



where we have used the initial condition 



i(k, 0)^,0)) = -%£(k + k'). 



Therefore we have 
A 



(10) 



(11) 



a(t) 



(27r) d 



d d ~ l k\ dfc||pr (2Tfe i +2fei) *~ 2fc » t ~ ( " o/3)fcib( * ) 



(27r) dV 2t y r[(d-l)/2] 



dk ± k d ± - 2 e-^ '*!+»!)*, (12) 



where r' = r + ^f^r-- Substituting (^ 7 ) 1 / 2 fc ± = x, this 
equation simplifies to 



(2ir) dy 2t J r[(d-l)/2] v ' 6 t 

dx x d-i e -(:S+* l )2tT>- 



(13) 



To solve this equation, we shall consider different dimen- 
sions separately. At d = 2, the saddle point is at the 
origin. Therefore we obtain 



a{tf = B {t + 



u b(t) 
6~T 



)/*. 



(14) 
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where Bq = 
we have 

2t 2 a(t) 



A 2 /4ir. Differentiating this equation once, 



da(t) uBq 



dt 
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a(t) + Bqt - 2ta(t) 2 



(15) 



The numerical solution of this equation in terms of three 
constants ci, C2 and C3 is 



a{t) = ci + + c 3 



logt 
i 



(16) 



In the large time limit, the dominant time dependence is 
of the form a(t) ~ t~ x l 2 . Finally, since the factor 2ir' 2 
in the argument of the exponential approaches infinity as 
t — » 00, the saddle point approximation turns out to be 
an exact one. 

Let us consider d = 3. Starting with (|l3|), we again re- 
place the integral by the maximum value of the integral. 
This occurs at 



As expected, the scaling of Ly remains unaltered from 
the Gaussian level. 

To summarize, we have studied the large n limit of a 
generalized n component driven diffusive model with an- 
nealed randomness at zero temperature. This model also 
corresponds to the order parameter conserved dynamics 
of uniaxial ferromagnets with dipolar interaction. This 
model is exactly soluble in the large n limit. The scaling 
of the length scale transverse to the drive is only modi- 
fied due to an additional logarithmic correction. This is 
quite different from the multiscaling behavior observed 
before in conserved dynamical models. 

I thank H. W. Diehl for the hospitality at Universitat 
Gesamthochschule Essen, where part of this work was 
performed. 



-2c ± y/4c 2 + 16c 
8c 



(17) 



where Cq = 2r' 2 t. As f -> 00, one can approximate 



1/2cq. Therefore in the large time limit, 



a(t) 



A 



4(2vr) 3 / 2 t 



exp[-(l/2 + l/4c )] 



(18) 



Therefore at d = 3, a(t) scales as 1/t at large time. To 
obtain the relevant length scale, we rewrite 4>i(k,t) as 



Mk,t) = i (fc,O)e- fc ' t e [( rt+ ( ll / 3! ) b ( t )) 2 / 4t l{ 1 -« fe / fc ^ 2+1 ) 2 > (19) [ 6 ; 



Apart from the length scale in the parallel direction grow- 
ing as iii ~ t 1 / 2 , we find that the length scale in the 
direction perpendicular to the driving force acquires a 
logarithmic correction in addition to the expected scal- 
ing behavior. At d = 3, the length scale in the transverse 
direction scales with time as 



L 



t l/4 



(l + ulnt/Qt 1 / 2 ) 1 / 2 



t 1 ^ for t 



(20) 



This is somewhat different from the multiscaling observed 
before, where there are two different scaling lengths. Here 
due to the anisotropy originally present in the system, 
the transverse and longitudinal length scales L± and Lu 
respectively, have different scaling properties with time 
even at the Gaussian level. L±, due to the non-Gaussian 
interaction, relevant only in the transverse direction, has 
modified scaling behavior due to a logarithmic correction 
in addition to the scaling expected at the Gaussian level. 



ir 

[8] 
[9] 

[10 

[11 

[12 
[13 
[14 

[15; 

[16 
[17 



B. Schmittmann and R. K. P. Zia in Phase transitions 
and critical phenomena edited by C. Domb and J. L. 
Lebowitz (Academic, London, 1995), Vol XVII 
T. Halpin-Healy and Y. C. Zhang, Phys. Rep. 254, 215 
(1995); S. Mukherji and S. M. Bhattacharjee, Current 
Science 77, 394 (1999) 

S. Katz, J. L. Lebowitz and H. Spohn, Phys. Rev. B 28, 
1655 (1983) 

M. Q. Zhang, Phys. Rev. A 35, 2266 (1987) 

S. Katz, J. L. Lebowitz and H. Spohn, J. Stat. Phys. 34, 

497 (1984) 

H. K. Janssen and B. Schmittmann Z. Phys. B 64, 503 
(1986) 

J. L. Valles and J. Marro, J. Stat. Phys. 49 89 (1987); J. 
Marro and J. L. Valles, J. Stat. Phys. 49, 121 (1987) 
B. Schmittmann, Euro. Phys. Lett. 24, 109 (1993) 

B. Schmittmann and R. K. P. Zia, Phys. Rev. Lett. 66, 
357 (1991) 

A. J. Bray, Adv. Phys. 43, 357 (1994) 

A. Coniglio and M. Zannetti, Euro. Phys. Lett. 10, 575 

(1989) 

C. L. Emmott and A. J. Bray, Phys. Rev. E 59, 213 
(1999) 

Humayun, Phys. Rev. Lett. 68, 



1559 



J. L. Cardy, J. Stat. Phys. 44, 567 
Siggia and H. A. Rose, Phys. Rev. A 



A. J. Bray and k. 
(1992) 

K. -t. Leung and 
(1986) 

P. C. Martin, E.D 
8 , 423 (1973) 

A. Aharony, Phys. Rev. B 8, 3363 (1973) 
J. Zinn- Justin, Quantum Field Theory and Critical Phe- 
nomena, International series ol monographs on Physics, 
Clarendon Press, Oxford, 3rd Edition, 1996 



4 



